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Asteroid Impact Risk Modeling 


“pats t=" Initial asteroid ¢ Fragment-Cloud Model (FCM) 
eae Sia propenles Models energy deposited in the 
atmosphere during entry and breakup 
e Energy deposition used to estimate 
airburst altitudes and ground damage 
e FCM results can also be matched to 
observed meteor light curves 
e Infer pre-entry asteroid properties 


e Investigate different breakup 
characteristics 


Entry & breakup 
energy deposition 


Airburst 


e Guide model refinements 
e Bound parameter ranges 
e Current effort is expanding FCM to 


represent varied initial asteroid 
structures 
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Fragment-Cloud Model (FCM) 
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Flight integration: 
dm/dt = -0.5p,;,V?Ao 
dv/dt = p,;.v’AC p/m — gsin@ 
d0/dt = (vV/(R;+h) — g/v)cosé 
dh/dt = vsin@ 


Fragmentation condition: 
Pairv’ > Strength (S) 


Fragment strengths increase 
with decreased size 
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Clouds broaden and slow 
under common bow shock 
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Cloud vs Fragment Contributions 


e Energy deposition computed as =. =o 
change in total kinetic energy of all 
fragment and cloud components as a 
function of altitude (kt/km) 

e FCM energy deposition mechanisms: 

e Debris clouds deposit the bulk of the 


Altitude (km) 


energy as they rapidly spread he: r 2 vo! 
: . Energy Deposition (kt/km) 
e Fragments serve to distribute the release 
of varied cloud masses es individual Cloud Generations 


e Large clouds released higher up produce 
broad, gradual flares 

e Small clouds released lower down 
produce sharper spikes 
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Observed Meteor Modeling 


e FCM approach combined with variable initial structures can match a 
range of energy deposition profiles and features from observed meteor 


light curves. 
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Chelyabinsk Modeling Approach 
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° Initial entry and mass estimates ¢ FCM parameters varied: 
(Popova et al. 2013): Initial rubble fragments/debris 
° 19.8 m diameter ¢ Initial aerodynamic strengths 
° 19.16 km/s at 18.3° ¢ Number and mass distribution of 
° 3.3 g/cm? meteorite-based density fragments per break 
¢ Vary FCM inputs to match energy e Cloud mass fraction per break 
deposition profile from observed light * Strength scaling exponent (a) 
curve (Brown et al., 2013). e Ablation coefficient (0) 
¢ Parameters are set as initial inputs and » Cloud dispersion coefficient (Cyisp) 


are not tuned along the entry. 


Energy Deposition from Light Curve 


Onsenved Right Sule (data from Brown et al, 201 3) 


FCM Breakup Modeling 
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e Cases using meteorite-based 
density (3.3 g/cc) consistently 
exceeded peak energy deposition 
by ~50% 

e Lower initial mass and bulk 
density better matched peak 
energy deposition 

¢ 9-10M kg (reduced from 13M kg) 

e Bulk density ~2.2—2.6 g/cc 

e Maintained 3.3 g/cc material density 

e Gives macroporosity ~21—33% 
compared to meteoritic density 


e Consistent with macroporosity 
ranges between fractured bodies 
(15-25%) or rubble piles (30-70%) 
(Britt et al., 2003, Asteroids III). 
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Pre-Entry Structure 
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_ 40 4 e Able to match all three 
™ main flare features using 
35 \ either a rubble pile, or a 


monolith with regolith. 
e Rubble pile provides 
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=_ = ¢ Regolith blow-off or 
< 35 rubble debris 
A ° Disruption ~0.5-0.6 MPa 
= 30 ° ~0.2-0.3% of initial mass 
<x e Not as well reproduced 
25 by shedding of small 
fragments 
20 ae 
15 
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e Stronger initial rubble: 
e ~2.5% of initial mass 
e Breakup at ~15.5 MPa 
e Successive 
fragmentation persisting 
below main flare 
e Relatively high strength 
scaling a (0.3-0.5) 
e 2-6 fragments per break 
with close (but not 
identical) mass splits 
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_ 40 & e Landed fragment 
= masses within ballpark 
35 \ of estimates from 


discovered falls 
e Total fallen fragment 
mass 5000-6500 kg 


25 =< * 4000-6000 kg estimated 
— in Popova et al. 2013 
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e Reducing cloud 
dispersion and/or 
ablation coefficients 
better matches 
width/slope of main peak 

¢ Dispersion Cys, = 1.5- 
2.5 (vs. 3.5 baseline) 

e Ablation o = 4-8e-9 kg/J 
(vs. 1e-8 baseline) 
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Summary & Future Work 
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e Developed FCM capability for modeling breakup and energy deposition of different 
asteroid structures 


e Produces realistic variety of energy deposition features, enabling very good matches to observed 
meteors. 


e Demonstrated how we can use those matches to make inferences about asteroid characteristics. 
e Found potential parameter refinements for modeling debris clouds. 


e Risk assessment applications 


e Analytic approach efficient enough to run the large numbers of cases needed for probabilistic risk 
assessments, yet variable enough to represent a wide range of potential asteroid structures. 


e Provides a way move beyond the typical point-source estimates and incorporate the different 
energy deposition rates into ground damage estimates. 
e Ongoing and future development: 
e Initialize rubble pile distributions using inverse power law distributions of sizes. 
e Refine parameters for cloud spread rates using hydrocode simulations 


e Explore effects of varied energy deposition profiles on ground damage compared to point- 
sources estimates (CFD simulations) 


e Automate a curve-matching optimizer to enable more comprehensive inference and thorough 
exploration of the parameter space 
e Paper available online: L.F. Wheeler et al., Icarus (2017), A fragment-cloud model for 
asteroid breakup and energy deposition. https://doi.org/10.1016/).icarus.2017.02.011 


May 18, 2017 L. Wheeler, PDC 2017 Page 14 


